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ABSTRACT 

Aims. Short phases of coeval powerful starburst and AGN activity during the lifetimes of the most massive galaxies are predicted by 
various models of galaxy formation and evolution. In spite of their recurrence and high luminosity, such events are rarely observed. 
Finding such systems, understanding their nature, and constraining their number density can provide key constraints to galaxy evolu- 
tionary models and insights into the interplay between starburst and AGN activities. 

Methods. We report the discovery of two sources at z=3.867 and z=3.427 that exhibit both powerful starburst and AGN activities. 
They benefit from multi-wavelength data from radio to X rays from the CFHTLS-Dl/SWIRE/XMDS surveys. Follow-up optical and 
near-infrared spectroscopy, and millimeter IRAM/MAMBO observations are also available. We performed a multi-wavelength analy- 
sis of their spectral energy distributions with the aim of understanding the origin of their emission and constraining their luminosities. 
A comparison with other composite systems at similar redshifts from the literature is also presented. 

Results. The AGN and starburst bolometric luminosities are ~10'^ Lq. The AGN emission dominates at X ray, optical, mid-infrared 
wavelengths, and probably also in the radio. The starburst emission dominates in the far-infrared. The estimated star formation rates 
range from 500 to 3000 Mo/yr. The AGN near-infrared and X ray emissions are heavily obscured in both sources with an estimated 
dust extinction Av>4, and Compton-thick gas column densities. The two sources are the most obscured and most luminous AGNs 
detected at millimeter wavelengths currently known. 

Conclusions. The sources presented in this work are heavily obscured QSOs, but their properties are not fully explained by the stan- 
dard AGN unification model. In one source, the ultraviolet and optical spectra suggest the presence of outflowing gas and shocks, and 
both sources show emission from hot dust, most likely in the vicinity of the nucleus. Evidence of moderate, AGN-driven radio activity 
is also found in both sources. Based on the estimated stellar and black hole masses, the two sources lie on the local Mbh - A^iu/ge 
relation. To remain on this relation as they evolve, their star formation rate has to decrease or stop. Our results support evolutionary 
models that invoke radio feedback such as the star formation quenching mechanism, and suggest that such a mechanism might play a 
major role also in powerful AGNs. 
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1. Introduction time le.g iHartwicS 12004 iHasinger et all l2005h and the cor- 

relation in nearby galaxies betwe e n black hole ( BH) m ass 
The similarity between the star formation history and the ^^^^^ jFerraresd l2002t iGebhardt et alJ IM in- 

space density of active galactic nuclei (AGNs) over cosmic ^-^^^^ jj^^j ^^^^ formation in a galaxy is related to the 



growth of its BH (but see [Shiel ds et all |2008[). The o retical 
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Although these models claim to b e successful iri reproducing 
and explaining several observations (Hopk ins et al.l2005btl2006t 
[Li et al. 2007; Chakrabarti et al. 2007 more observations are 
necessary to test their predictions, constrain their parameters, 
and provide a ph ysical base to some of their assumptions (see 
e.g. lMarullietan i2008). 

Previous studies of SMG populations and high-z starbursts 
have revealed the presence o f obscured AGN activity in ~30- 
40% of these systems (e. g . iPage et all 12 004; Chapm an et al 



2005 



Stevens et aP '2005^, 'Alexan der et al.. ,2005g; ,Yan et al 



2007 : IPo'pe et 31.120 08). giving support to a hnk between the two 



activities. However, it is not clear how star formation and BH 
growth are linked, what their relative timescales are, and whether 
they influence each other directly through a feedback mecha- 
nism or indirectly, e.g. by consuming the available cold gas. The 
number of well studied composite systems is still low, and it is 
often difficult to separate and quantify the contribution from dif- 
ferent energy sources (see e.g. the variet y of models proposed to 
expla i n the properties of F 10214-H472 4; Rowan-Robinson et al.l 
[1991 IrbphtzeTanilOOet fef stathioul l2Cm In order to quan- 
tify their contribution to the bolometric luminosity, the multi- 
wavelength spectral energy distributions (SEDs) of such sys- 
tems need to be measured and modeled. Because of the difficulty 
of obtaining a full multi-wavelength coverage, AGN, starburst 
galaxies, and composite systems with well sampled SEDs and 
known spectra are often used to derive correlation between bolo- 
metric luminosities and measurements at specific wavelengths o r 
of spectral features (e.g. IS aiina e"taL[|2008l: [Polletta et aLliioOSl) . 
In case of powerful starburst galaxies, like SMGs, the AGN is 
often highly obscured and requires either ultra deep X-ray ob- 
servatio ns (A lexander et aI 2005a), spectropolarimetric obser- 
vations ( [Goo drich et alj^ 1 996h. or mi d-infrared (MIR) spectro- 
scopic observations (P ope et al.ll2008l) to be revealed. 

In order to investigate the hnk between star formation and 
AGN activity, it is thus important to study systems where both 
activities are taking place, and collect multi-wavelength mea- 
surements to constrain the starburst, i.e. the star formation rate 
(SFR), and AGN, i.e. accretion rate, luminosities, as well as their 
Eddington ratio, and stellar and gas masses. 

Here, we investigate the properties of two rare millimeter 
(mm) bright obscured QSOs at high-z discovered in a wide 
multi-wavelength survey. Spectroscopic data as well as the full 
SED, from X ray to radio wavelengths are available for both 
objects and are analyzed to constrain the origin of their lumi- 
nosity and investigate their nature. Throughout this paper, we 
adopt a flat c osmology with Hp - 71 kms"' Mpc Qm=0.27 
and Qa=0.73 dSpergel et al.ll2003h . 



2. Target selection and observations 

In this work, we analyze the properties of two sources that 
were selected as g drop-out sources, thus z-4- candidate, in 
the 0.9deg2 CFHTLS survey Dl[3 (RA = 02''26'", Dec = 
-04°30') and as very bright MIR sources in the XMM-LSS 
field of the Spitze r Wide-area InfraR ed Extragalactic Legacy 
survey (SWIR]^; Lonsdale^etafl [20031). The sources lAU of- 
ficial names are SWIRE2 J022550.67-042 142.2 (SW022550 
hereinafter), and SWIRE2 J022513.90-043419.9 (SW022513 
hereinafter). We applied a modified version o f the d rop-out se- 
lection technique developed by Steidel et"aL[ d 19991) for z >3- 
4 Lyman-break sources. In addition to the optical ugr colors 



and the source extension, this selection technique takes into ac- 
count the 3.6yum-4.5yum color (Bergeron et al., in prep.; see 
also Siana et al. 2008). No other sources with these properties 
were found in the 0.9 deg^ CFHTLS-Dl/SWIRE field. Because 
of their large 24/im fluxes (F24;,m > 1 mJy), exceptionally red 
infrared (IR) SEDs (aju <-2.5 where Fy oc v''"' over the 
observed wavelength range 3-24 yum), and high-z, we selected 
both sources for observations at 1.2 mm with the Max Plan ck 
Milhmeter Bolometer (MAMBO) array dKrevsa etalJll998l) at 
the Institut de Radioastronomie Millimetrique (IRAM) 30m tele- 
scope. Follow up spectroscopic observations were carried out 
and broad-band photometric data from X ray to radio wave- 
lengths were already available in the field. 



2.1. Optical and infrared imaging 

Broad-band photometric data at optical and IR wavelengths 
are available for both sources from various surveys. Optical 
data in 5 broad-bands, ugriz, were provided by the CFHTLS 
survey Dl (data release TOOO-©. Near-IR (NIR) data in the 
J and K bands were obtained from the UKIRT Infrared 
Deep Sky Survey (U KIDSSB data release 3; iDve et all [20061; 
Lawrence et aLluOOTh . Infrared data in all seven Spitzer bands, 
from IRAC (3.6, 4.5, 5.8, and 8.0yum; Fazio et al. 2004) and 
MIPS (24, 70, and 160jum; IRieke et al.1 12004 ). are available 
from the SWIRE survey. The SWIRE data correspond to the 
latest internal catalogs that will be released as part of the 
Data Rel ease 5 (DR5; for details on the data reduction see 
ISurace e t al. 2005). The total measured magnitudes and fluxes 
in each band are reported in Table [T] 



2.2. MAMBO observations 

Observations at the IRAM 30m telescope were carried out dur- 
ing March- April and November 2006, using the 117 element 
version of the MAMBO array operating at a wavelength of 
1.2 mm (250 GHz). We used the standard on-off photometry ob- 
serving mode, chopping between the target and sky at 2 Hz, and 
nodding the telescope every 10 or 20 s. On-off observations were 
typically obtained in blocks of 6 scans of 16 or 20 lOs-subscans 
each, and repeated in later observing nights. The atmospheric 
transmission was intermediate with T(1.2mm)=0.1-0.4. The ab- 
solute flux calibration was established by observations of Mars 
and Uranus, resulting in a flux calibration uncertainty of about 
20%. 

On average, the noise of the channel used for point-source 
observations was about 35-40 mJy/ Vf/beam, where t is the 
exposure time in seconds, consistent with the MAMBO time 
estimator for winter conditions. The two sources were ob- 
served as part of a larger program targeting more than 100 
starbur st, AGN, and compos ite sources selected in the SWIRE 
fields dLonsdale et al.l [20081 Polletta et al., in prep., Fiolet et 
al., in prep.). The resulting 1.2 mm fluxes are 4.7+0.8 mJy and 
5.5±0.7mJy for SW022550, and SW022513, respectively. 

2.3. Optical spectroscopy ofSW022550 

An optical spectroscopic observ ation of SW02255 was carried 
out with the AAOmega system dSharp et al.ll2006h on the 3.9-m 
Anglo- Australian Telescope, on September 27, 2006, as part of 



' http : //www. cfht. Hawaii .edu/Science/CFHTLS/ 

^ |http : // Swire . ipac . caltech . edu/ swire/swire . html] 



http://terapix . iap . f r/rubrique . php?id-rubrique=24Tl 



http : //www . ukidss . org/ 
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a multi-object spectroscopic program of QSOs in the CFHTLS 
(Petitjean et al., in prep., Yu et al. in prep.). The 580V and 385R 
grating were used in the blue and red arms of the spectrograph, 
respectively, yielding R ~ 1300. The total integration time was 
3h, divided into six 30min exposures. 

Data reduction was performed using the AAOmega data re- 
duction pipeline software DRCONTROL. The two dimensional 
images were flatfielded, and the spectrum was extracted (using 
a gaussian profile extraction), wavelength calibrated and com- 
bined within DRCONTROL. The final spectrum was flux cali- 
brated using the broad-band photometric measurements in the g 
and r bands after deconvolving the spectrum with the filter band- 
passes. This calibration assumes that the source optical emission 
has not significantly varied between the spectroscopic and pho- 
tometric observations. 

The optical spectrum of SW022550 shows several strong 
emission lines, e.g. Lya, N V /11240, C IV A1549, the O VI 
/11035 doublet (see Figure [TJ. The continuum emission is weak 
(<2cr), and a drop in the continuum emission is visible at rest- 
frame wavelengths below 1216A due to intergalactic medium 
(IGM) absorption. The bottom panel of Figure [T] shows the 
decrease in the continuum emission at rest-frame wavelengths 
/1<1216A. Absorption features due to the Lyman or forest are 
also visible. Since the line spread function (LSF) derived from 
the lamp spectra is consistent with a Gaussian, we used such 
a profile to fit the emission lines. From the fits we measured 
the total line flux, the central wavelength, the FWHM in A, and 
the flux of the local linear continuum. The FWHM in kms ' 
and the rest-frame equivalent width (Wa) in A are then derived 
from these parameters after correcting for instrumental resolu- 
tion. The instrumental FWHM is 3 A at /1<5500A, and 5 A at 
/l>5500 A. The line flux, central wavelength, corresponding red- 
shift, rest-frame FWHM and rest-frame W,i of the main visible 
lines are reported in Table |2] Uncertainties are determined us- 
ing a Montecarlo method and assuming a Gaussian noise given 
by the r.m.s. associated with the continuum where the fit is per- 
formed. A broad component, in addition to a narrow one, is mea- 
sured in correspondence of the Lya, N V A124Q, Si IV A1391, 
and C IV /11549 lines. All these lines show an asymmetric com- 
ponent, similar to a blue wing. To fit such an asymmetric profile, 
in some cases, e.g. Lye, andN V /11240, we fit multiple Gaussian 
components. Based on the central wavelength of the optical lines 
listed in Table|2] we estimate a redshifts of 3.867±0.009. 

The emission lines are dominated by their narrow 
(FWHM- 1000 kms ') components wit h FWHMs typ i cal of 
type 2 QSOs (< 1500 kms"'; see e.g. [Baldwin etal] Il988t 
iNorman et al.l 120021; iMainieri et al.l |2005'). Broad components 
are also present with FWHM~2000- 9000 kms ' as observ ed in 
type 1 QSOs (>2000 kms see e.g. lRichards et al.ll2004 . The 
equivalent widths of the emission lines in the optical spectrum 
of SW022550 are larger than observed in type 1 QSOs, and 
more similar to those observed in other type 2 QSOs (see e.g . 
SMM 02399-0136 CDFS-263 or CXO 52; llvison et al.l[T99i 
IMainieri et al.ll200llStern et al.ll2 002). In Figure^ we compare 
the optical spectrum of S W022550 with a composite spectrum of 
a large sample of type 1 QSOs and with the spectra of four type 
2 QSOs at high-z from the literature. We show the spectra nor- 
malized at the Lya peak on the top panel, and at the C IV /II 549 
peak in the bottom panel. The composite type 1 QSO spectrum 
corresponds to the median composite spectrum of 2204 QSO 
spectra from SDSS (Vanden Berk et al. 2001). The type 2 QSOs 
have been selected because of the availability of their optical 
spectra from the literature and high-z. They are CDFS-202, and 
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Fig. 1. AAOmega optical spectrum of SW022550. The expected 
emission lines are labeled. The horizontal dashed line corre- 
sponds to a null flux. The bottom panel shows the same spectrum 
but smoothed with a kernel of 7 pixels to better illustrate the con- 
tinuum emission. The top dotted grey line represents the mean 
continuum level at rest-frame /l=1275-1380 A, and the bottom 
soHd grey curve the mean continuum level at /l=1050-l 170 A. 



CDFS-263 from CDFS (ISzokolv et alJ |2004 . an d SW104406 
and SW104409 from the SWIRE/Chandra survey (P oUetta et aP 
20061 This comparison shows the wide range of line ratios that 
can be observed in type 2 QSOs. The main differences reside 
in the strength of the N V /1 1240, and C IV /1 1549 emission 
lines, and the blue wings observed in the main emission lines. 
We will discuss below two possible explanations for the strength 
of the N V /11240, and C IV /11549 lines, i.e. high metalicities 
and shocks. Interestingly, there is another object in the type 2 
QSO sample from the literature with similar asymmetric lines, 
i.e. SW104409. This source is also characterized by a similar 
optical-lR SED to SW022550, with an unusual optical blue con- 
tinuum with strong AGN emission lines, and red optical-MlR 
colors. A possible explanation for the presence of blue wings is 
extinction due to dust mixed with infalling or outflowing ionized 
gas (Osterbrock 1989). Since the SED in both cases indicate that 
the AGN is obscured, it is also plausible that the blue contin- 
uu m and the broad blu e wings are due to scattered light (see e.g. 
Zak amska et al.ll2006l) . 

SW022550 is characterized by large N V .11240/C IV .11549 
and N V AUAQ/Ue II /11640 flux ratios. These flux ratios ai-e 
often used to estimate the gas metalicity in QSOs, from ei- 
ther the broad or the narrow line region (BLR and NLR; see 
Nagao et al. 2006, and references therein). High value s usu- 
ally imply high metalicities (see e.g. iHamann & Ferlandlll993l; 
IVernet et al.l l2001l 'Norman et al.l j20q2l). Ba sed on ±s photoion- 
ization model in Hamann & Ferland ( 1993b . we derive a metalic- 
ity Z-4Zq. However, not all the emission line ratios of the spec- 
trum of S W 022550 are explained by NLR pho toionization mod- 
els (see e.g. iNagao et al1l2006t iGroves et alll2004 . A possible 
cause of this discrepancy is the presence of broad components. 
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Table 1. Multi-wavelength data of SW022550 and SW022513 
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Fig. 2. Top panel: Optical spectra normalized at the flux at 
1216 A in the rest-frame of SW022550 (black thick Hne), of 
the type 2 QSOs CDFS-263 (c y an lin e) and CDFS-202 (ma- 
genta line) from Szokoly et akl (l2004. and SW104406 (blue 
Une), and SW104409 (red hne) from iPolletta etalj (|2006|) . and 
medium composite spectrum of SDSS type 1 QSOs (green line; 
[Vanden Berk et al. 2001). Bottom panel: As top panel, but nor- 
maUzed at the flux at 1549 A rest-frame. The main emission lines 
are labeled and marked with dotted vertical lines. The horizon- 
tal dashed line corresponds to a null flux. All spectra have been 
smoothed for clarity. 



It is also possible that the emission lines are produced by shocks 
in addition or rather than by photoionization. This hypothesis is 
supported by the strength of some lines, like C IV /1 1549, O VI 
/11035, and N V AllAQ. For example, C IV /11549 is predomi- 
nantly produced in the cooling zone of shocks and it is expected 
to be stronger in shocks than in photoionization processes. 
Indeed, the C IV ^1549/He II ^1640, N V ^1240/He II /II 640, 
and N V /I1240/C IV /1 1549 flux ratios like those observed in 
SW022550, are higher than predicted by photoionization, and 
can be explained by shocks dAUen et al.l 2008), or by the co- 
existence of shocks and AGN photoionization , without requir- 
ing la rgely super-solar metalicities ( Moy & Rocca-Volmerangj 
I2OO2I) . Also the presence of O VI /(1035 can be considered as 
evidence of shocks, although other explanations are also pos- 
sible. To confirm whether SW022550 is characterized by high 
metalicity and better constrain its value, or whether shocks are in 
part responsible for the observed emission lines, it would be nec- 
essary to measure other emission lines, especially at rest-frame 
optical wavelengths. 

In summary, the optical spectrum of SW022550 shows spec- 
tral features that are also observed in other obscured (type 2) 
AGNs at high- z, but are not expl ained by the standard AGN uni- 
fication model d Antonuccil 1 9931) . In particular, a high metalicity 
or the presence of shocks and a significant scattered fraction at 
rest-frame UV wavelengths are necessary to fully explain the ob- 
served spectrum. Note that these hypothesis can not been fully 
tested with the available data and other scenarios are not ruled 
out. 



2.4. Near-infrared spectroscopy 

Both targets were observe d with the ISAAC instru- 
ment (iMoorwood et al.l 1 19981) on Antu (VLT-UTl) in low- 
resolution (LR) mode, using the 1024x1024 Hawaii Rockwell 
array of the Short Wavelength arm. Observations were carried 



mag or Flux^ 


Proj ect/Tele scope 


SW022550 


SW022513 


Fo.5-2/SfV 


XMDS/XMM 


<2.5 


1.1+0.8 


F2-10/SfV 


XMDS/XMM 


<20 


7.8+4.1 


m„ 


CFHTLS/CFHT 


25.211+0.077 


>26.00 


m^' 


CFHTLS/CFHT 


22.944+0.007 


24.432+0.033 


m^' 


CFHTLS/CFHT 


21.648+0.003 


23.287+0.014 


m,v 


CFHTLS/CFHT 


22.949+0.004 


22.851+0.012 


m. 


CFHTLS/CFHT 


22.408+0.021 


22.503+0.027 


p; 


UKIDSS/UKIRT 


3.6+0.4 


3.1+0.6 


F^ 


UKIDSS/UKIRT 


7.7+0.6 


12.7+0.6 




SWIRE/Spitzer 


15±I 


16+1 




SmRE/Spitzer 


23+1 


20+1 


Fs.S^m 


SWIRE/Spitzer 


<58 


28+4 


Fg.O^m 


SWIRE/Spitzer 


239+7 


180+5 




SWIRE/Spitzer 


3.32+0.02 


2.35+0.02 


F^o^lm 


SWIRE/Spitzer 


<24 


<24 




SWIRE/Spitzer 


<I26 


<126 


Fl.2/)?/j? 


IRAM/MAMBO 


4.70+0.77 


5.53+0.72 


F20rai 


VIMOS/VLA 


0.14+0.02 


0.35+0.02 


FsOcni 


VIMOS/GMRT 


0.41+0.04 


0.62+0.04 



" X ray fluxes are in 10 ergs cm "s '. Optical magnitudes are total 
magnitudes in the AB systems in the ugriz filters from the CFHTLS 
survey Dl. UKIDSS J and K fluxes, and IRAC (3.6-8.9yum) fluxes are 
in p]y. MIPS (24-160 pm), MAMBO ( 1 .2 mm) and radio (20 cm) fluxes 
are in mjy. X-ray upper limits correspond to 3cr, optical upper limits 
correspond to 80% completeness, IR upper limits correspond to 5cr. 
The r adio data are from the VIMOS VLA survey CBondi et al..2003l 
l2007h . 



out with a 1" slit with two grating blaze angles, approximately 
covering the SH {A ~ 1.42-1. 82//m), and SK (A ~1.82-2.5yum) 
bands. The nominal resolutions were 500, and 450 for the two 
cases, respectively. The observations of the two targets were 
caiTied out in visitor mode on September 13'''-14''', 2007. 
Total integration times were 8 1 min and 66 min in the SH band, 
and 84 min and 42 min in the SK band on SW022550, and 
SW022513, respectively. 

Each observation was splitted in sets of 180s exposures, 
adopting a nodding and jittering pattern for optimal background 
subtraction. The nod throw was set to 60-80 ", in order to avoid 
conflict with the reference star; the jitter width was 10". 

Data reduction was carried out in the standard manner, us- 
ing the IRAF environment, including ABBA sky subtraction. 
Wavelength calibration was based on OH sky lines. Atmospheric 
extinction correction and flux calibration were obtained with 
standard telluric stars, carefully selected from the ESO database 
in order to have known NIR magnitudes. For this purpose, sev- 
eral B-type telluric standards were observed each night (before 
and after a target was observed) at air masses within 0.1-0.2 
of the air mass of the target observations. The rebinned ISAAC 
SK-band spectra of the two sources are shown in Figure [3] and 
discussed below. 

2.4.1 . ISAAC SK Spectrum of SW022550 

Neither continuum nor lines are detected in the SH-band spec- 
trum of SW022550. Two emission lines are detected in the SK- 
band spectrum at /i23560A, and 24060A yielding a z -3.85 if 



' The package IRAF is distributed by the National Optical 
Astronomy Observatory which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agree- 
ment with the National Science Foundation. 
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Fig. 3. The top panels show the ISAAC SK spectra of SW022550 (left panel) and SW022513 (right panel). The grey curves 
correspond to the spectrum smoothed with a kernel of 2 pixels, and the black curves to the smoothed spectrum with a kernel of 
10 pixels. The bottom panels show the SK band sky spectrum. The location of the main expected lines are labeled. 



interpreted as H/3, and [O III] /14959, respectively. The contin- 
uum in the SK-band is not detected and it is estimated to be 
<3.8xlO"'^ ergs cm"^ s"' A"' (corresponding to 3cr). The iden- 
tification of these lines is tentative due to the high sky level at 
these wavelengths and the uncertain wavelength calibration at 
A >2Apm. The [O III] /i4959 line shows a blue tail probably 
due to a poor sky subtraction. At the expected location of the 
[O III] A5QQ1 line, the spectrum is negative, suggesting that the 
sky subtraction might have been too high at these wavelengths. 
The H/3/Lya flux ratio is 0.03+0.1, which is consistent with the 
range of values observed in type 2 AGNs, i.e. from 0.014 to 
0.03, while type 1 AGNs are usually characterized by higher val- 
ues (0.07-0.125). Note that because of the poor sky subtraction 
where the lines are observed (see bottom panel of Figure [3]l, it is 
possible that a broad Hfi component is also present. Assuming 
a FWHM of 5000 kms ' for a broad H/S we estimate an up- 
per limit to the flux < 16x10"'^ ergs cm"^ s"'. This value is sig- 
nificantly higher that the predicted flux based on the measured 
Lya flux and the H/3/Lya flux ratio observed in type 1 AGNs. 
Thus, we cannot rule out the presence of a broad H/S line in 
the spectrum of SW022550. The [O III] A4959fH/3 flux ratio 
is 6.3±2.1. This is also closer to the ratios observed in type 2 
AGNs, -2.3, than to those observed in type 1 AGNs, -0.04. The 
upper hmit to flie li/3 flux yields a [O III] /14959/H/3 flux ratio 
>0.56, still higher than in type 1 QSOs. In summary, in spite 
of its low signal-to-noise, we conclude that the (rest-frame) op- 
tical spectrum of SW022550 is more similar to the spectra of 
type 2 AGNs, than to those of type 1 AGNs. The ISAAC spec- 
trum shows another interesting property. The H/3, and [O III] 
/14959 lines are slightly blueshifted with respect to the optical 
emission lines. Assuming that such a difference is real and that 
the source is at z=3.867 as derived from the optical spectrum, 
the blueshifted lines could be explained by outflowi ng gas with 
a velocity ~510 kms"' (see similar cases in ISaiina et al. 2008). 
Suc h a velocit y is consistent with those assumed in shocks mod- 
els dAUen et aL.200&) . 



2.4.2. ISAAC SK Spectrum of SW022513 

Neither continuum nor lines are detected in the SH-band spec- 
trum of SW022513. In the SK-band spectrum, two emission 
lines at /lo/, 5 =2 1954 A, and 22 160 A are well detected, while 
the continuum is not detected. We thus estimate an upper limit 
to the continuum <2. 03x10""* ergs cm"^ s"' A"' (correspond- 
ing to 3cr), and to the emission lines equivalent widths. The two 
lines are interpreted as [O III] /1/14959,5007, constraining the 
redshift to be z=3.427±0.001. A weak emission feature is ob- 
served at the expected wavelength corresponding to the H/3 line. 
Since the line is affected by a bright sky line which falls exactly 
at its expected location, and it is narrower than the spectral reso- 
lution, we cannot constrain it. The [O III] /14959/Hy8 and [O III] 
/15007/Hj6 flux ratios are, respectively, >0.97 and >2.48. These 
values are much higher than those observed in type 1 AGNs 
(i.e. 0.04 and 0.15) and more similar to those observed in type 2 
AGNs (2.3 and 5.5) and expecte d from photoionization models 
of the narrow line region (NLR; lOsterbrockll 1 9891; iGroves et al.l 
2004). Thus, it is quite likely that such a feature is the Hj3 line. 
In summary, the ISAAC spectrum of SW022513 shows narrow 
emission lines with flux ratios consistent with those of type 2 
AGNs. 

The lines fluxes, Wa, and FWHMs of both spectra are re- 
ported in Table|2]The SH-band observations cover the rest-frame 
wavelength range -2900-3700 A for SW022550, and -3200- 
4 100 A for SW022513. At these wavelengths, the only feature 
that might be detected is the [O II] /i3727 emission line. The 
[O II] /13727 Une is typically associated with star formation ac- 
tivity ([Gallagher et al.ll 19891) . but also with AGN activity. Since 
our sources are not detected in the SH-band, and an atmospheric 
feature is present at the expected [O II] A3727 location in the 
two sources, we cannot set any constraints on their [O II] A3121 
emission. 

From the spectra we estimate that the upper limits to the 
MR continuum in SW022513 and SW022550 are <33juJy, and 
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<61 yuJy. These values are about 3 and 8 times higher than the 
K-band measured fluxes. Thus, the spectroscopic measurements 
are consistent with the broad band photometric data. 

2.5. Radio data 

Both targets are detecte d in the VIMOS VLA survey at 
1.4 GHz, and at 6 10 MHz (" Bondi et al.ll2003HCiliegi etani2005t 
iBondi et al. 2007). The measured fluxes are reported in Table[T] 
In the VLA 1.4 GHz observations, SW022513 appears extended 
with a major axis of 4.1" and a minor axis of 1.8", correspond- 
ing to a projected line ar size of 4.2kpcx 2. 1 kpc, and a position 
angle (PA) of 161.3° ( Bondi et al.ll2003l) . In the other radio ob- 
servations the sources are unresolved. Assuming a power-law 
model Fy oc V'' for the radio emission, the radio spectral indeces 
ar are -1.3 for SW022550, and -0.7 for SW022513. Because of 
its radio spectral index, SW022550 can be co nsidered an ultra- 
steep spectrum source (USS, i.e. Or < - 1 : Roettgeri ng et alj 
[l994, and references therein). Such a steep spectrum implies an 
AGN origin of the observed emission. In case of SW022513, the 
spectral index is consistent with what is observed in star forming 
galaxies and radio-quiet AGNs dCiliegi et al.l2003l) . however the 
extended emission and the large radio flux strongly suggest also 
an AGN origin of its radio emission. A more detailed discussion 
of the radio properties of our sources is presented in § 13.31 

2.6. X ray data 

The two targets are also covered by 20 ks observations 
with XMM-New ton from the X MM- Newton Me dium Deep 
Survey (XMDS; Chiappetti et aTll2005t iTaier et al.i r2007). Only 
EPIC/MOS data are available as both sources fall in a bad col- 
umn of the EPIC/pn detector SW022513 is at 4.5' from the aim- 
point and marginally detected with 29±12 total counts (12±8 
at 0.5-2 keV, and 17±9 at 2-10 keV). The hai'dness ratio is 
HR=0.18^y J, which corresponds to an effective gas column den- 
sity Nh-(1^o5)x10^'* cm"^ assuming a power-law model with 
photon index r=1.7, and Galactic and int rinsic photo-electric 
absorption with Nh'^°'= 2.61x10^" cm"^ jPickev & LockmanI 
11990 ). More d etails on the X ray dat a and derived quantitie s 
can be found in lChiappetti et al.1 (l2005h . and lTaier et al] (l2007h . 
SW022550 is at >11' from the aim-point and is not detected, 
thus we can only give an upper limit to its X-ray flux. The X-ray 
fluxes of both sources are reported in Table[T] 

Based on these X ray measurements, the estimated broad 
and hard X ray luminosities of SW022513 are ~6xl0"^^ ergs s"', 
and 5x10'*'* ergs s"', respectively. The absorption-corrected 
broad and hard X ray luminosities are 11x10'*'* ergss"', and 
6x10'*'* ergss"'. For SW022550, we estimate an upper limit 
to the absorbed luminosities in the broad and hard bands of 
2.0x10'*^ ergs s"', and 1.8x10'*^ ergs s"', respectively. Based on 
the estimated X ray luminosity of SW0221513 and the relation- 
ship between the [O III] /t5007 and the hard X ray lumi nosity in 
type 2 AGNs dMulchaev et alJll994tlNetzer et alj|2006l) . we ex- 
pect an [O III] A5001 luminosity of 3-9x10'*^ ergss"*, instead 
of 2x10'*'* ergs s"' as observed (see TableO. Such a discrepancy 
might be due to the non validity of the X ray/[0 III] /15007 re- 
lationship for our sources, to an underestimation by a factor of 
100 of the X ray luminosity, or to the contamination from star 
formation to the [O III] A5QQ1 flux. Probably, all these factors 
play a role in explaining this discrepancy, but the only one that 
could explain such a large difference is the underestimation of 
the intrinsic X ray luminosity. A higher X ray luminosity would 



also be required to obtain a MIR/X-ray luminosity ratio con- 
sistent with those observed in other AGNs (see § 13. It . Based 
on these considerations, it is quite probable that SW022513 is 
a Compton-thick QSO. Since SW022550 is even fainter in the 
X-rays, and more luminous in [O III] /15007 emission and in 
MIR luminosity, it is also quite plausible that SW022550 is a 
Compton-thick QSO. 

The full radio-Xray SEDs of the two targets are shown in 
Figure|5] and all the available fluxes are listed in TabIe[T] 

3. Multi-wavelength spectral energy distributions 

The SEDs of SW022550, and SW022513 are shown in Figuresi] 
and |5] Both sources are characterized by large MIR fluxes 
(F24^m-2-3 mJy), show red power-laws in the NIR (l-5jum in 
the rest-frame; ajR<-2.5, where Fyoc v"'"), and large IR/optical 
flux ratios (Log(vF24^m/yFz)-l.5). Al l these properties are 
signatures typical of obscured QSOs jWeedman et al.l 120061 : 
iPoUetta et al.]|2008l) . 

3.1. Comparison witli otiier liigti-z QSOs 

We compare the UV-mm SEDs of SW022550, and SW022513 
with those of well studied mm-detected non radio-loud AGNs at 
high-z from the literature in Figure |4] The vast majority of mm- 
detected QSO at high-z are unobscured or type 1 QSOs with 
MIR SE Ds consistent with the median template of unobscured 
QSO (e.g.'Hi nes et al.l2006l) . We represent th eir SEDs in panel a 
of FigurelHwith the median QSO template by lElvis et 10(11994 . 
The majority of these type 1 QSOs with strong mm fluxes, i.e. 
>3mJy, as our sources, show a significant excess in the far- 
IR (FIR) with respect to the optical emission predicted b y the 
median QSO template (IWang et al.ll2007l: iHao et alJl2008D . The 
SEDs of 4 QSOs with a FIR excess are also shown in panel a 
of Figure m Our sources show similar mm luminosities, but 10- 
100 times lower optical fluxes, consistent with being obscured. 
In Figure |4] we also compare our sources SEDs with those of 
3 mm-detected obscured QSOs with available MIR data, (1) the 
lensed system F 10214-1-4724, a Seyfert 1 galaxy at z=2.86, but 
that i s also heavily absorbed i n the X rays (panel b; llvison et"an 
II998 UAIexanderetal.ll2005bh: (21 SMM 02 399-0136. a type 2 
QSO atz=2.803 (panel c: llvison et al.ll998h : and (3) CXO GWS 
JI4174I.9-H52 2823 (CXOJ1417 hereinafter) at z=1.15 (panel e; 
iLe Floc'h et a l. 2007). With the only exception of CXOJ1417 
that is the least luminous sources of the sample, all the remain- 
ing cases are characterized by mm luminosities consistent with 
those of our sources. All the obscured AGNs shown in panels 
b-d in Figure |4] show flat NIR SEDs with similar luminosities. 
Since the maximum of stellar light in vFy and the minimum of 
AGN light are at NIR wavelengths (I-1.6/zm in the rest-frame; 
ISawickil l2002t ISanders et 111119891) . it is at these wavelengths 
that we can expect the maximum contribution to the optical-IR 
SED from the host galaxy. The similarity of the observed NIR 
SEDs might thus be explained by a significant contribution from 
stellar emission. All SEDs are characterized by red MIR col- 
ors, consistent with hot dust thermal emission, but such a com- 
ponent is more luminous and redder in our sources than in the 
sources from the literature. Our sources are also systematically 
more luminous at optical wavelengths, and this is mainly due to 
the strong emission lines. In summary, our targets have SEDs 
that are intermediate between those of type 1 QSOs (panel a), 
and those of obscured QSOs (panels b-d). In absence of high 
spatial resolution data and MIR spectroscopy, we cannot deter- 
mine the origin of the observed differences. Assuming that the 
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Line 


^obs 


z 




Flux 


FWHM" 


Comments 




(A) 




(A) (10 


ergs cm"^ s"') 


(kms-') 




SW022550 


Ly/? 


4997.4±0.8 


3.876 


41±3 


1.35±0.11 


935±180 


Narrow 


[O VI]i^l031,1038 


5018.2±0.6 


3.867 


68±3 


2.25±0.10 


777±81 






5042.9±0.8 


3.858 


152±4 


5.03±0.14 


1820+114 




Lyff 


5928±0.1 


3.875 


322±5 


13.52±0.19 


1041 + 13 


Narrow 


N V^1240 


6022±0.7 


3.856 


311±9 


13.19±0.39 


2531+99 


Broad 




6044±0.7 


3.874 


146±5 


6.19±0.20 


1354+74 


Narrow 


Lyo- 


5920±1.8 


3.868 


164±15 


6.88±0.65 


4060+157 


Broad 




5985±2.2 




496±23 


20.98±0.96 


8968+343 


Broad 


Si IV /II 397 


6778±4.3 


3.852 


833±37 


7.38±0.33 


4581+318 






6820±2.4 


3.861 


281 ±22 


2.49±0.20 


2375+172 




N IV] Ame 


7242±2.4 


3.873 


41±4 


1.95±0.18 


1987+239 




CIVA1549 


7527±1.5 


3.859 


328±8 


16.13±0.38 


3690+101 


Broad 




7542±0.8 


3.869 


91±4 


4.47±0.18 


1235+88 


Narrow 


He 11^1640 


7970±15 


3.860 


84±8 


2.75±0.26 


3098+768 






23560 


3.847 


>37* 


1.4±0.5 


<510" 


ISAAC 


[O III] ^4959 


24060 


3.852 


>239'' 


9.1±0.4 


<498'' 


ISAAC 


SW022513 




21532+5 


3.429 


>120'' 


<2.5 


<55r 


ISAAC 


[O III] ^4959 


21954±2 


3.427 


>297'' 


6.1±0.5 


288+58 


ISAAC 


[O III] .15007 


22160±1 


3.426 


>829'' 


16.9±0.6 


660+34 


ISAAC 



All measurements are based on single Gaussian fits to the emission lines assuming a flat (in F^) continuum. "Corrected for instrumental resolution. 
*Since the continuum is not detected we estimate an upper limit to the equivalent width, W^, assuming a 3cr upper limit to the continuum. 
The continuum is <3.78xl0-'** ergscm^^ s"' A"' for SW022550, and <2.03xlO-"* ergscm"' s"' A"' for SW022513. ''The FWHM is fixed to the 
spectral resolution because it cannot be correctly determined due to the low signal-to-noise. 




Fig. 4. Rest-frame UV-mm SEDs of the selected targets (black full circles), SW022550 (top panels), and SW022513 (bottom pan- 
els) , compared with the SEDs of other mm-detected AGNs, (a) a med ian type 1 QSO te mplate normalized at the mm luminos- 
ity jElvis et al .' 1994), and 4 mm-detected type 1 QSOs with FIR excess dWang et al.i20 07, purple symbols), J033829.31+002156.3 
(purple plus signs), J075618.14-H410408.6 (purple tr iangles), J092721.82-H200123 7 (purple diamond s), and J104845.05+463718.3 
(purple cr osses), (b) F 10214 -H4724 (green crosses: iRowan-Robinson et al.lll993l; iTeplitz et al.ll2006l) . (c) SMM 02399-0136 (red 
diamonds: iMson et al.|[T998l) . and (d) CXOJ1417 (brown triangles: iLe Floc'h et al.ll2007h . 



intrinsic AGN SED is the same in all sources, the redder MIR 
SED and the higher optical flux imply a higher AGN luminos- 
ity and more extinction in our sources than in the mm-detected 
obscured QSOs from the literature. 

In Figure |5] we compare our sources radio-X ray SEDs 
with those of high-z AGNs with similar multi-wavelength cov- 



erage. We include the mm-detected sources CDFS-263, and 
CXOJ1417 because they are also detected at X ray wavelengths. 
CXOJ1417 is also detected at radio wavelengths. We also in- 
clude C DF-202, a X ray dete cted heavily absorbed type 2 QSO 
at z=3.7 (iNorman et al.l(2002l) . and Mrk231, a BAL QSO which 
is heavily absorbed in the X rays and is hosted by a pow- 
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Rest— frame Log(i/) (Hz) Rest— frame Log(i/) (Hz) 

Fig. 5. Rest-frame radio-X ray spectral energy distributions of the selected targets (black full circles), SW022550 (left panel), and SW022513 
(right panel) compared with the SEDs of X ray detected heavily absorbed AGNs from the literature, CDF-202 at z=3.7 (magenta circles; 
Gorman etal. 2 002; Miller etal. 2008), CDFS-263 at z=3.66 (blue squares; Mainieri et al. 2005), CXOJ1417 at z=1.15 (brown triangles; 
iLe Floc'h et al.ll2007 '). and Mrk231 at z=0.042 (purple solid line). The rectangles represent the soft and hard X ray fluxes and associated un- 
certainties. Downward arrows represent Scr upper limits. 



Rest-frame Logfi/) (Hz) Rest-frame Loq(u) (Hz) 

15 14 13 12 15 14 13 12 




-10 12 3 -10 12 3 

Rest-frame Log(X) (^im) Rest-frame Log(X) {/J.m) 



Fig. 6. Rest-frame UV-mm spectral energy distributions of the selected targets (cyan full circles), SW022550 (left panel), and 
SW02251 3 (right panel). The red dotted curves represent the starburst models normalized at t he observed 1 .2 mm flux from the 
library in lSiebenmorgen & Kriigell ( l2007l) . The dashed blue curves represent torus models from lHonig et alj (12006 ) normalized at 
24jum. The green solid line is a type 1 AGN template normalized at the observed J-band (1.25yum) flux. The solid purple curves 
represent the total of the starburst, torus, and type 1 AGN templates which yield the 50 hestx^, where the;\f^ is computed using 7 
data points from the J-band to the 24 jum, and do not overpredict the observed fluxes and upper limits at FIR wavelengths (see text 
for more details). The models are used to estimate the luminosities reported in Table[3] 



erful starburst galaxy. We do not consider F 10214+4724 be- 
cause its^£ra^_S2ectruri^ likely dominated by starburst emis- 
sion ([Alexander et al.l l2005b). The absorption-coiTected hard X 
ray luminosities of CDFS-263, CXOJ1417, and CDF-202 are 
7.6xl0'*4^ 2.4xl0'*4^ 3.3x10^4 ergss"', respectively. The X 
ray luminosity of SW022513 is ~6xl0'*'* ergss"'. Unless, the 
X ray luminosity of SW022513 is largely underestimated (see 
§ 12. 6t , the cuiTent estimate is consistent with those measured in 
the literature sources. In addition to being more luminous at all 
wavelengths, our sources show a redder MIR emission than the 
CDFS sources. Interestingly, such a red MIR emission is also ob- 
served in other heavily obscured QSOs at high-z which are also 
characterized by high er MIR/X ray lumin osity ratios (~13^0) 
than classical AGNs (Polletta et al.ll2008h . For classical AGNs 
the MIR/X ray luminosity ratios (vLg^m/La-iofcv) range from 



about 0.3 to 13 dPoUetta et al.ll2007l l2008l) . but for SW022513 
and CXOJ1417 tiiis ratio is about 60. Such a prominent MIR 
component might imply either that dust heating is more efficient, 
that the dust covering factor is particularly high, or that the in- 
trinsic X-ray luminosity is underestimated. 

3.2. SED modeling and star formation rate 

Our targets have optical and MIR properties typical of obscured 
QSOs, but their large mm fluxes suggest that they are also ex- 
periencing intense star forming activity. Indeed, emission from 
AGN-heated circumnuclear dust is expected to peak around 10- 
30yu m and no t to contribute significantly in the FIR (see e.g. 
Gra nato et al.ll2 004; Honig et al. 2006; Fritz et al. 2006). A pow- 
erful starburst provides a more likely explanation, but such a hy- 
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pothesis cannot be confirmed because of the lack of detections 
at rest-frame wavelengths A=6-250fim, contrarily to the cases 
for F 10214+4724 and CXOJ1417. Although current data do not 
constrain the origin of the mm emission, AGN model predictions 
and the similarity with F 10214+4724 and CXOJ1417 favor a 
starburst origin. 

Assuming that the mm flux is due to a starburst component 
and that the MIR fiux is mainly due to the AGN, we model the IR 
(2-1200 ;um) observed SED of our targets by combining a star- 
burst and a torus model. To model the starburst component, we 
use a library of 7208 starburst models (ISiebenmorgen & Kriigell 
12007 ). To model the AGN component, we use a set of 44 
clump y torus models with various inclinations from lHonig et al.l 
( |2006|) . These models were used to fit the MIR SEDs and spec- 
tra of a sample of MIR-luminous obscured QSOs with avail- 
able MIR spectra (Polletta et al. 2008). The starburst compo- 
nent is normalized at the observed 1.2 mm flux, and the torus 
models at the observed 24 //m flux. The sum of the starburst 
and torus models is then compared to the SED of our targets. 
Since none of the models is able to reproduce the optical data, 
and since the optical data indicate that there is a contribution 
from strong emission lines, typical of AGNs, an additional com- 
ponent is added to the starburst+torus models represented by 
an AGN template. For this component we use a typ e 1 QSO 
template (Hatziminao glou et all 120051: iPolletta et all 1200 7) be- 
cause it reproduces well the observed optical broad-band data. 
The type 1 QSO template is normalized at the J-band (1.25 jum) 
observed flux. Note that such a normalization corresponds to 
~l.4% of what would be obtained if we normalized the type 
1 QSO template at the observed 24 /im flux. Such a small frac- 
tion implies that the AGN optical light might be suppressed by 
more than 70%, which corresponds to a foreground extinction of 
Ay =4.6. Since this component is not reddened, its origin is more 
likely due to a small (1.4%) fraction of nuclear light that escapes 
without suffering obscuration, or to nuclear scattered light (see 
e.g. Zakamska et al. 2006). 

After combining all the starbursts models, with all the torus 
models and with the type 1 QSO template, we then reject all 
combinations that overpredict the fluxes at 70, and 160/im. This 
criterion rejects about 40% of all models. We then select 50 mod- 
els with the lowest x^, computed by comparing the predicted and 
observed fluxes in 7 bands, from 1 .25 /vm to 24 yum. The best 50 
models are shown in Fig. |6] The choice of 50 models yields a 
wide enough range of solutions to characterize all acceptable 
models and the luminosities produced by each component. A 
higher number of models yields either poor fits or does not in- 
crease significantly the range of SEDs already represented by the 
50 best models. 

The preferred starburst models contain an old stellar popula- 
tion component, correspond to starbursts of intermediate sizes, 
i.e. 3 kpc, and show a wide range of ratios, ~40-90%, of the OB 
stars luminosity with respect to the total luminosity. The pre- 
ferred torus model for both sources, defined as the most often 
chosen among the best 50 models, corresponds to a torus with 
axis at only 15° from the line of sight with an additional cold 
dusty absorber along the line of sight. Such a model explains 
the prominent MIR emission as radiation from hot dust in the 
torus inner walls, and the red MIR colors as extinguished torus 
radiation by a cold absorber (for details on the AGN models see 
iPoUetta et al.]|2008l) . 

From the set of acceptable models we derive the AGN and 
the starburst bolometric (0. 1-1000 ;um) luminosities, and the IR 
(8-1000 ;um) starburst luminosity. The derived luminosities for 
all acceptable models are reported in Table |3] 



Table 3. SW022550 and SW022513 luminosities 



Parameter" 



Method* 



SW022550 SW022513 



Log(L(IR)) 

Log(Li.4GH3) 

Log(L(H)) 

Log(M.) 

Log(L(6yum)) 



SK07 Starburst 
PL 

SK07 Starburst 
L(H) (S07) 
H06 Torus 
SK06 Starburst 
H06 Torus 



12.5- 13.3 
25.50 

11.2-11.9 
11.2-11.9 
13.1-13.2 

12.6- 13.4 
13.4 



12.5-13.2 

25.38 
11.3-11.6 

11.3- 11.6 
12.9-12.9 

12.4- 13.2 
13.1 



" L(IR) is the 8-1000 yum luminosity in Lq derived by integrating 
the model under Method. L140H; in the monochromatic radio lumi- 
nosity at 1.4 GHz in the rest-frame in W Hz '. L(H) is the monochro- 
matic luminosity in Lq at 1.6/jm in the rest-frame derived from the 
starburst models. M, is the stellar mass derived from L(H) assuming 
th e average L(H)/M. rat io derived from the high-z radio galaxy sample 
in lSevmour et al.l ( 120071) . L(6//m) is the rest-frame luminosity at 6/jm 
in Lq derived from the torus models. L^^ is the 0. 1-1000 /jm luminos- 



T AGN 



is the 0. 1-1000 yum 



ity derived from the starburst models in Lq. 
luminosity derived from the torus models in Lq. 

* The SK07 starburst mod els are from the library 

in ISieb enmorgen & Kriigell (l2007l) . The H06 torus models are 
from Hon ig et al. U2006I) . PL stands for power-law model with spectral 

index a,-, Fy oc v"% where ar=-l.3 for SW022550, and -0.7 for 
SW022513. 



Assuming that the estimated IR luminosity is produced 
by a starburst, the derived SFRs, g iven by SFR(MQ/yr) = 
L(IR)/(5.8xlO''L0) (iKennicuttI Il998h. are -500-3000 Mo/yr 
for both sources. The estimated SFRs are consistent with 
those estimated in other radio-quiet obsc ured QSOs at high- 
z (e.gllvison e tal. 1998; Efstathiou 2006: IStevens et al.1 120051; 
iPoIIetta et alT2008i;iSaiina et al.ii2007b,) . 

3.3. The origin of the radio luminosity 

Among the various feedback mechanisms that have been 
suggested as star formation regulators in high-z galax- 
ies, the radio feedback mode is probably the most pop- 
ular and successful in reproduc i ng numerous obse r vables 
through simulations (ISilk & ReesI 1 19981 ICroton et al.1 120061; 
ICattaneo & Tevssie^ l2007h . Because of the predicted short 
timescale associated with the radio activity and of the difficulty 
of identifying AGN-driven radio activity in non radio-loud ob- 
jects, there is only limited evidence of radio feedback at play on 
star-formation in AGNs, especially at high-z, and at high lumi- 
nosities. The best supporting observations are the detection of 
AGN-driven outflows in high-z radio-galaxies (Nesvadb a et al.l 
2006, ,2007b. a), and of moderate AGN-driv en radio activity in 
hig h-z 5piteer- selected starburst and AGNs (ISaiina et al.ll2007at 
PoIIettal I2OO8I) . Finding further evidence of AGN-driven ra- 
dio activity, especially in massive high-z star forming galax- 
ies, would provide additional support for the radio feedback as 
quenching mechanism. Here, we analyze the origin of the radio 
emission in our two sources. 

The measured radio fluxes of SW022550, and 
SW022513 correspond to 1.4 GHz rest-frame luminosi- 
ties of 10^5 50 w Hz- 'and lO^^ ^8 WHz-', respectively, or 
SFRs~6OOOM0/yr, if powered by star formation. These high 
luminosities would imply the presence of AGN-driven radio 
emission in the local universe (an AGN is classified radio-loud if 
PiAGHz > 10-^^ WHz \ corresponding to a SFR ^lOOMo/yr; 
Condon 1992). However, at high-z, intense starburst episodes 
with SFRs of thousands of Mo/yr can be expected. It is thus not 
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straightforward to infer an AGN origin of the radio emission 
in sources at high-z based on the estimated radio luminosities. 
Although a broad multi-wavelength coverage is available for 
both sources, some traditional diagnostic methods would not be 
reliable in probing an AGN origin of the radio emission. The 
commonly diagnostic based on large radio/optical flux ratios 
cannot be applied because our sources suffer from heavy optical 
obscuration. Another diagnostic is provided by the radio/infrared 
flux ratio, which is based on the well established far-IR/radio 
correlation (ICondonlll992h . The average FIR-to-radio emission 
ratio q = log(L(FIR)/3 .75 x 1 ' ^ W) -log(Lv. , agh .I^ Hz" ' ) 
is eq ual to 2.34 in local star forming galaxies Vim et al.l 
l200l . where L(FIR)=L4o-i20/.«„ and L(F1R)/L(1R)=0.56±0.1 
for the starburst models used in this work. Slightly lower 
values have been fou nd in sources at higher z, q-2.Ql in z ~2 
SMGs dKovacs et all [2006). an d t7=2.21±0.02 in radio-quiet 
MIR selected z -1-2.5 AGNs ( Saiina et al.l|200 8'). As the AGN 
contribution to the radio emission increases, q decreases. For 
example, in radio loud AGNs q is typically close to zero or 
negative. The q values derived for our sources using the FIR 
luminosities obtained from the starburst models (see Table O 
are ~0.8-1.6 for both sources. These low values favor an AGN 
origin for the observed radio emission. This is also supported 
by the steep radio spectral index in SW022550 and the extended 
radio size of SW022513. 

In summary, we find some evidence of an AGN origin for the 
radio emission of our two sources. The associated radio activity 
is moderate, as it is only slightly in excess of what is typically 
observed in radio quiet quasars and much less than in radio loud 
quasars. 

3.4. AGN and starburst contribution to tiie bolometric 
luminosity 

The SEDs of our sources indicate that their MIR emission is 
dominated by AGN-heated dust. Using the best torus models de- 
scribed in § l3.2l and shown in Fig.|6] we estimate a lower limit to 
the AGN bolometric luminosity by integrating the torus model 
in the 0. 1-1000 yum wavelength range. The derived AGN bolo- 
metric luminosities, also reported in Table [3] are lO'^ '^L©, and 
10'^ 1 Lo for SW022550 and SW022513, respectively. To derive 
the starburst bolometric luminosity, we do the same integration 
but using only the starburst models. Based on these approxima- 
tions, the AGN contribution to the bolometric luminosities is 
-50-87% in SW022550, and 43-83% in SW022513. The broad 
range of values is due to the uncertainty on the starburst IR lu- 
minosity. The AGN contribution to the system bolometric lumi- 
nosity is thus at least 40%, and it can be more than 80%. 

3.5. Host stellar luminosity 

The stellar emission, in Fy, from a host galaxy p eaks in the NIR , 
typically at 1 .6 fim in the rest-frame (or H-band) jSawic kil2002h . 
The ratio between NIR luminosity and stellar mass is character- 
ized by little dispersion, especially when similar star formation 
histories are assumed. The estimated stellar masses can be 30% 
lower in case of younger stellar populations. Thus, we can as- 
sume that the NIR luminosity is a proxy of the stellar mass or 
luminosity of the host galaxy. Since in our objects the AGN 
optical and NIR emissions are obscured, we can directly esti- 
mate the host NIR luminosity from the observed SED. Indeed, 
the best-fit SED models show that the NIR emission is domi- 
nated by stellar light. We estimate the host, and thus stellar, H- 



band host luminosities using the acceptable starburst templates 
described above. From the H-band luminosity, we estimate the 
stellar mass using the average L(H)/M, ratio derived for a sam- 
ple of high-z radio galaxies by ISevmour et al.l (l2007l) . Log(M,) 
= Log(L(H))-0.1±0.1, where M, is the stellar mass in Mq and 
L(H) is the monochromatic H-band lumin osity in Lq. Note that 
this relationship assumes a Kroupa] (1200 Ih initial mass function 
(IMF). The estimated stellar masses are 1.3-6.5x10" M© for 
SW022550, and 1.5-2.9x10" Mq for SW022513. The estimated 
H-band luminosities and stellar masses are reported in Table [3] 
Because of the uncertainty in separating the host and the AGN 
contributions to the H-band luminosity, the stellar masses should 
be considered as upper limits to the true stellar masses. The es- 
timated masses are of the sam e order as those of the most mas- 
sive systems at high-z (see e.g.'Sevmour e t al.]|2007l: Ifierta et al.l 
2007), and are thus among the most massive objects at z~3-4. 

4. Discussion 

4.1. Dust obscuration: geometry and distribution 

The two sources analyzed in this work exhibit properties usu- 
ally observed in type 2 AGNs. Their narrow line dominated 
spectra, low optical/IR flux ratios and faint X ray fluxes im- 
ply that they are obscured, but both sources show a blue opti- 
cal continuum which i s unexpected by the standard AGN unifi- 
cation model dAntonuc ci 1993), but not unusual for obscured 
AGNs (see e.g. SW104409; PoUetta et all |2006|) . Moreover, 
SW022550 shows some broad components in its optical spec- 
trum. In addition, both sources also exhibit a prominent MIR 
emission suggesting that the hottest dust component, likely at 
the dust sublimation radius, is in part visible (see preferred torus 
models in § 13. 2t . These hybrid properties suggest that a frac- 
tion of the nuclear light is visible. Since the SEDs of our sources 
cannot be simply reproduced by extinguishing the emission of a 
type 1 QSO, a more complex model is required to explain their 
properties. A proper modeling, that takes into account various 
geometries, dust properties, the intensity and the spectrum and 
luminosity of the heating source, and a radiative transfer treat- 
ment would be necessary to fully explain their hybrid properties. 
However, the uncertainty on determining the spectrum of the 
AGN component due mainly to the contribution from the host 
galaxy, and the lack of a MIR spectrum, would yield degenera- 
cies in the models. Thus, we limit our discussion to a qualitative 
comparison with some of the most popular models of obscur- 
ing dust in AGNs from the literature. A possible scenario is a 
clumpy dust distribution which produces a partial c overing of 
the nucleus ( Nenkova et"aLll2002l; iHonig et alJ l2006V This sce- 
nario is quite unlikely if the optical radiation comes from the 
nuclear region, and it is thus much smaller than the dust sub- 
limation radius. An alternative to an unobstructed line of sight 
to the nucleus is having a scattering medium, perhaps associ- 
ated with the NLR. The scattering medium can produce the blue 
optical light, and broad emission lines. If dust is mixed with 
the scattering medium, this might be responsible for the large 
MIR excess, while the nuclear light is obscured by the dusty 
torus jSchweitzer et alj |2008). Spectropolarimetric observations 
to measure the scattered component are planned for SW022550 
to investigate the latter scenario. 

4.2. SMG-QSO evolutionary link 

Our targets show simultaneous powerful starburst and AGN ac- 
tivity. They are thus good candidates for transitioning objects be- 
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tween the starburst and the QSO stage predicted by the most pop- 
ular ev olutionary models (" Sanders et al.lll988t [Alexander et alj 
I2005d) . Here, we consider their properties to make predictions 
on their evolution. 

From the AGN bolometric luminosity (see Table O, and as- 
suming that the AGN emission is Eddington limited (Lt,,! - 
LEdd), we derive a lower limit to the BH masses of SW022550 
and SW022513 of ^ 7.4 x lO** Mq and 3.7x10'* Mg, respec- 
tively. Sincejn^st^ftheqim accrete below their Eddington 
limit dMcLure & Dunlopl2004l) . it would be more realistic to as- 
sume a lower Eddington ratio, but this would imply even higher 
BH masses. Note that these estimates are not quite accurate 
since we are estimating the AGN bolometric luminosity from 
the reprocessed thermal emission which is highl y dependent on 
the p roperties of the reprocessing dust (see e.g. ' Marconi et al.l 
12004 ). Assuming an accretion efficiency of 10%, the derived ac- 
cretion rates are -17 Mq yr"' for SW022550 and 8.5 Mq yr"' 
forSW022513. 

Assuming the estimated BH masses and the stellar masses 
derived in § 13.51 we find that these sources lie on the local M^h- 
^hidge relationship (e.g. Marconi & Hunt 2003). Thus they are 
more similar to other obscured QSOs at high-z, than to SMGs 
with AGN activity or type 1 AG Ns which are both offset from the 
local MBH-Mbuhe relatio nship ([Alexander et all 120081; iPoUettal 
120081: ICoppin et aLll200 8'). In order to keep these objects on the 
local MsH-^buige relation as they evolve, the on-going star for- 
mation process and BH growth will have to stop, or continue 
with analogous rates (SFRsM). Since accretion rates of thou- 
sands, or even hundreds of M^/yi are not sustainable, the latter 
scenario can occur only if the SFR decreases significantly. In ei- 
ther cases, we predict a quick decrease in the SFR of our sources. 
The presence of shocks and outflowing gas and the moderate 
AGN-driven radio activity might be the expression of an AGN- 
induced feedback that will cause the predicted SFR to decrease. 
Such a mechanism i s predicted in numerous current evolutionary 
models ( Silk 20P5i ICroton et al.1 .2006), but it is more often in- 
voked in AGNs at lower redshifts and luminosities. Our results 
suggest that ra dio feed back might be also important in high-z 
QSOs (see also lSajina et al. 2008; Polletta 2008) . 

5. Summary 

We analyze the multi-wavelength spectral energy distributions 
of two obscured QSOs at high-z discovered in the CFHTLS- 
Dl/SWIRE survey, SWIRE4 J022550.67-042 142.37 (z=3.867) 
and SWIRE4 J022513.92-043420.24 (z=3.427). Both sources 
benefit from multi-wavelength data available in the field and 
from mm data from MAMBO and optical and NIR spectroscopy 
obtained as part of a dedicated follow up program. 

Their large mm fluxes, >4mJy, imply IR luminosities 
^^q12.5-i3.3Lq, and SFRs~500-3000Mo/yr The two sources 
also show powerful AGN activity, mostly seen in the MIR 
(L(MIR)aga' >1O'^L0) and in the optical emission lines, and 
probably also at radio wavelengths. Their MIR emission is 
mainly due to AGN-heated hot dust and corresponds to a lu- 
minosity comparable to or higher than the starburst luminosity. 
The MIR SEDs are redder and more luminous than observed in 
other mm-detected obscured QSOs (see § 13. H and Figure|4]i. 

Their optical and NIR spectra are characterized by faint or 
absent continuum emission and strong emission lines typical of 
AGNs. The optical spectrum of SW022550 suggests either a 
high metalicity, Z ~ 4-Zq, or the presence of shocks. The compar- 
ison between the rest-frame ultraviolet (UV) and optical spectra 
indicates the presence of outflowing material with a velocity of 



-500 kms~^ The red and faint NIR emission compared to the 
MIR emission indicates that the AGN is extincted by dust with 
Av>4. The faint X ray emission, the large hardness ratio, and 
the small X ray/[0 III] /15007, and X ray/MIR luminosity ratios 
suggest that their X ray emission might be absorbed by a gas col- 
umn density >10^'* cm"^, and thus that they are Compton-thick 
QSOs. 

The sources show some hybrid properties in between those 
of type 1 and type 2 AGNs, i.e. a broad component in some 
UV emission lines, blue UV colors, and a significant emission 
from hot dust, likely at the sublimation radius. We suggest that 
these properties are explained by the presence of scattered light 
at UV wavelengths, and that a hot dust component might be vis- 
ible, perhaps associated with the NLR and with the scattering 
medium. 

The scattering scenario will be tested with planned spec- 
tropolarimetric observations on the VLT of SW022550. With 
those data we will investigate whether the UV emission of 
S W022550 is dominated by a scattered component, and measure 
the scattered fraction. 

The composite (starburst and AGN) nature of our sources 
suggests that they might be experiencing the transition phase be- 
tween the starburst a nd the QSO stage pred icted by some evolu- 
tionary models (e.g. lDi Matteo et al]|2005h . Since both sources 
already lie on the local Mbh - Mbuige relation, we predict a 
quick decrease in their SFR in order to remain on such a relation 
as they evolve. Both sources are characterized by AGN-driven 
moderate radio activity, providing support to evolutionary mod- 
els invoking radio feedback a s star formation quenchin g mecha- 
nism even in powerful QSOs (ISilkll2005l; ICrotonll2006l) . 

Acknowledgements. We thank the referee, R. Antonucci, for a very careful read- 
ing and helpful comments and suggestions. MP is grateful to S. Andreon for 
providing a list of clusters around our sources. MP is also grateful to M. Bondi 
for providing the VLA image of SW022513. MP thanks N. Fiolet for the fi- 
nal validation of the MAMBO measurements. MP is particularly thankful to 
N. Miller and C. Norman for providing the radio data for CDF-202, and for 
helpful discussions. MP thanks S. Honig for the permission to use the torus 
models. CSS is grateful to S. Croom for providing the AUTOZ code. MP ac- 
knowledges financial support from the Marie-Curie Fellowship grant MEIF-CT- 
2007-042111. The AAOmega observations have been funded by the Optical 
Infrared Coordination network (OPTICON), a major international collabora- 
tion supported by the Research Infrastructures Programme of the European 
Commission's Sixth Framework Programme. This work is based on observations 
made with the Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology under NASA contract 1407. We 
are thankful to the IRAM staff for their support with the observations and the 
data reduction. 

References 

Alexander, D. M., Bauer, F E., Chapman, S. C, et al. 2005a, ApJ, 632, 736 
Alexander, D. M., Brandt, W. N., Small, I., et al. 2008, AJ, 135, 1968 
Alexander, D. M., Chartas, G., Bauer, F E., et al. 2005b, MNRAS, 357, L16 
Alexander, D. M., Small, I., Bauer, F E., et al. 2005c, Nature, 434, 738 
Allen, M. G., Groves, B. A., Dopita, M. A., Sutherland, R. S., & Kewley, L. J. 

2008, ApJS, 178, 20 
Antonucci, R. 1993, ARA&A, 31, 473 

Baldwin, J. A., McMahon, R., Hazard, C, & Williams, R. E. 1988, ApJ, 327, 
103 

Berta, S., Lonsdale, C. J., Polletta, M., et al. 2007, A&A, 476, 151 
Bondi, M., Ciliegi, P, Venturi, T, et al. 2007, A&A, 463, 519 
Bondi, M., Ciliegi, P, Zamorani, G., et al. 2003, A&A, 403, 857 
Cattaneo, A. & Teyssier, R. 2007, MNRAS, 376, 1547 
Chakrabarti, S., Cox, T. J., Hemquist, L., et al. 2007, ApJ, 658, 840 
Chapman, S. C, Blain, A. W., Small, I., & Ivison, R. J. 2005, ApJ, 622, 772 
Chiappetti, L., Tajer, M., Trinchieri, G., et al. 2005, A&A, 439, 413 
Ciliegi, P, Zamorani, G., Bondi, M., et al. 2005, A&A, 441, 879 
Ciliegi, P, Zamorani, G., Hasinger, G., et al. 2003, A&A, 398, 901 
Condon, J. J. 1992, ARA&A, 30, 575 

Coppin, K. E. K., Swinbank, A. M., Neri, R., et al. 2008, MNRAS, 389, 45 



12 



Polletta et al.: Obscured and powerful AGN and starburst at z~3.5 



Croton, D. J. 2006, MNRAS, 369, 1808 

Croton, D. J., Springel, V., White, S. D. M., et al. 2006, MNRAS, 365, 1 1 
Di Matteo, T., Springel, V., & Hemquist, L. 2005, Nature, 433, 604 
Dickey, J. M. & Lockman, F. J. 1990, ARA&A, 28, 215 
Dye, S., Warren, S. J., Hainbly, N. C, et al. 2006, MNRAS, 372, 1227 
Efstathiou. A. 2006. MNRAS, 371, L70 

Elvis, M., Wilkes, B. J., McDowell, J. C, et al. 1994, ApJS, 95, 1 
Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 
Ferrarese, L. 2002, ApJ, 578. 90 

Fritz, J., Franceschini, A., & Hatziminaoglou, E. 2006, MNRAS, 366, 767 
Gallagher, J. S., Hunter, D. A., & Bushouse, H. 1989, AJ, 97, 700 
Gebhardt, K., Bender, R., Bower, G., et al. 2000, ApJ, 539, L13 
Goodrich, R. W., Miller, J. S., Martel, A., et al. 1996, ApJ, 456, L9 
Granato, G. L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 2004, ApJ, 
600, 580 

Granato, G. L., Silva, L., Monaco, R. et al. 2001, MNRAS, 324, 757 
Groves, B. A., Dopita, M. A., & Sutherland, R. S. 2004, ApJS, 153, 75 
Hamann, F & Ferland, G. 1993, ApJ, 418, 11 

Hao, C.-N., Xia, X.-Y, Shu-DeMao. Deng, Z.-G., & Wu, H. 2008, Chinese 

Journal of Astronomy and Astrophysics, 8,12 
Hartwick, F D. A. 2004, ApJ, 603, 108 

Hasinger, G., Miyaji, T., & Schmidt, M. 2005, A&A, 441, 417 
Hatziminaoglou, E., Perez-Fournon, I., Polletta, M., et al. 2005, AJ, 129, 1 198 
Hines, D. C., Krause, O., Rieke, G. H., et al. 2006, ApJ, 641, L85 
Honig, S. R, Beckert, T., Ohnaka, K., & Weigelt, G. 2006, A&A. 452, 459 
Hopkins, P R, Hernquist, L., Cox, T. J., et al. 2005a, ApJ, 630, 705 
Hopkins, P R, Hernquist, L., Cox, T. J., et al. 2005b, ApJ, 630, 716 
Hopkins, P R, Hemquist, L., Cox, T. J., et al. 2006, ApJ, 639, 700 
Ivison, R. J., Small, 1., Le Borgne, J.-F, et al. 1998, MNRAS, 298, 583 
Kennicutt, Jr, R. C. 1998, ARA&A, 36, 189 

Kovacs, A., Chapman, S. C, Dowell, C. D., et al. 2006, ApJ, 650, 592 
Kreysa, E., Gemuend, H.-R, Gromke, J., et al. 1998, in Proc. SPIE Vol. 3357, 

p. 319-325, Advanced Technology MMW, Radio, and Terahertz Telescopes, 

Thomas G. Phillips; Ed., ed. T. G. Phillips, 319-325 
Kroupa, P 2001, MNRAS, 322, 231 

Lawrence, A., Warren, S. J., Almaini, O., et al. 2007, MNRAS, 379, 1599 
Le Floc'h, E., Willmer, C. N. A., Noeske, K., et al. 2007. ApJ, 660, L65 
Li, Y., Hernquist, L., Robertson, B., et al. 2007, ApJ, 665, 187 
Lonsdale, C. J., Polletta, M., Omont, A., et al. 2008, ApJ accepted 
|arXiv:0810.1293| 

Lonsdale, C. J., Smith, H. E., Rowan-Robinson, M., et al. 2003, PASP 115, 897 
Mainieri, V., lUgopoulou, D., Lehmaim, I., et al. 2005, MNRAS, 356, 1571 
Marconi, A. & Hunt, L. K. 2003, ApJ, 589, L21 
Marconi, A., Risaliti, G., Gilh, R., et al. 2004, MNRAS, 351, 169 
Marulh, F, Bonoli, S., Branchini, E., Moscardini, L., & Springel, V. 2008, 
MNRAS, 257 

McLure, R. J. & Dunlop, J. S. 2004, MNRAS, 352, 1390 

Miller, N. A., Fomalont, E. B., Kellermann, K. I., et al. 2008, ApJS accepted 
|arXiv:0804.2375| 

Moorwood, A., Cuby, J.-G., Biereichel, P., et al. 1998, The Messenger 94, 7 

Moy, E. & Rocca-Vohnerange, B. 2002, A&A, 383, 46 

Mulchaey, J. S., Koratkar, A., Ward, M. J., et al. 1994, ApJ, 436, 586 

Nagao, T, Maiolino, R., & Marconi, A. 2006, A&A, 447, 863 

Nenkova, M., Ivezic, Z., & Elitzur, M. 2002, ApJ, 570, L9 

Nesvadba, N. P. H., Lehnert, M. D., De Breuck, C, Gilbert, A., & van Breugel, 

W. 2007a, A&A, 475, 145 
Nesvadba, N. R H., Lehnert, M. D., Eisenhauer, R, et al. 2006, ApJ, 650, 693 
Nesvadba, N. R H., Lehnert, M. D., Genzel, R., et al. 2007b, ApJ, 657, 725 
Netzer, H., Mainieri, V., Rosati, P, & Trakhtenbrot, B. 2006, A&A, 453. 525 
Norman, C, Hasinger, G., Giacconi, R., et al. 2002, ApJ, 571, 218 
Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and active galac- 
tic nuclei (Research supported by the University of California, John Simon 
Guggenheim Memorial Roundation, University of Minnesota, et al. Mill 
Valley, CA, University Science Books, 1989, 422 p.) 
Page. M. J., Stevens. J. A., Ivison, R. J., & Carrera, F J. 2004, ApJ, 611, L85 
Polletta, M. 2008, A&A, 480, L41 

Polletta, M., Tajer, M., Maraschi, L., et al. 2007, ApJ, 663, 81 
Polletta, M., Weedman, D., Honig, S., et al. 2008, ApJ, 675, 960 
Polletta, M., Wilkes, B. J., Siana, B., et al. 2006, ApJ, 642, 673 
Pope, A., Chary, R.-R., Alexander, D. M., et al. 2008, ApJ, 675, 1171 
Richards, G. T, Nichol, R. C, Gray, A. G., et al. 2004, ApJS, 155, 257 
Rieke, G. H., Young, E. T., Engelbracht, C. W., et al. 2004, ApJS, 154, 25 
Roettgering, H. J. A., Lacy, M., Miley, G. K., Chambers, K. C, & Saunders, R. 

1994, A&AS, 108, 79 
Rowan-Robinson, M., Efstathiou, A., Lawrence, A., et al. 1993, MNRAS, 261, 

513 

Sajina, A., Yan, L., Armus, L., et al. 2007a, ApJ, 664, 713 
Sajina, A., Yan, L., Lacy, M., & Huynh, M. 2007b, ApJ, 667, L17 



Sajina, A., Yan, L., Lutz, D., et al. 2008, ApJ, 683, 659 

Sanders, D. B., Phinney, E. S., Neugebauer, G., Soifer, B. T., & Matthews, K. 
1989, ApJ, 347, 29 

Sanders, D. B., Soifer, B. T., Ehas. J. H., Neugebauer, G., & Matthews. K. 1988, 

ApJ, 328, L35 
Sawicki, M. 2002, AJ, 124, 3050 

Schweitzer, M., Groves, B., Netzer, H., et al. 2008, ApJ, 679, 101 
Seymour, N., Stern, D., De Breuck, C, et al. 2007, ApJS, 171, 353 
Sharp, R., Saunders, W., Sinith, G., et al. 2006, in Presented at the Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 6269, 
Ground-based and Airborne Instrumentation for Astronomy. Edited by 
McLean, Ian S.; lye, Masanori. Proceedings of the SPIE. Volume 6269, pp. 
62690G (2006). 

Shields, J. C, Walcher, C. J., Boker, T, et al. 2008, ApJ, 682, 104 

Siana, B., Polletta, M. d. C, Smith, H. E., et al. 2008, ApJ, 675, 49 

Siebenmorgen. R. & Kriigel, E. 2007, A&A, 461, 445 

Silk, J. 2005, MNRAS, 364, 1337 

Silk, J. & Rees, M. J. 1998, A&A, 331, LI 

Spergel. D. N., Verde, L., Peiris, H. V, et al. 2003, ApJS, 148, 175 
Springel, V, White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629 
Steidel, C. C, Adelberger, K. L., Giavalisco, M., Dickinson, M., & Pettini, M. 

1999, ApJ, 519, 1 
Stem, D., Moran, E. C, Coil, A. L., et al. 2002, ApJ, 568, 71 
Stevens, J. A., Page, M. J., Ivison, R. J., et al. 2005, MNRAS, 360, 610 
Surace, J. A., Shupe, D. L., Fang, F., et al. 2005, Spitzer Science Center 

California Institute of Technology 
Szokoly, G. P, Bergeron, J., Hasinger, G., et al. 2004, ApJS, 155, 271 
Tajer, M., Polletta, M., Chiappetti, L., et al. 2007, A&A, 467, 73 
Teplitz, H. I., Ai-mus, L., Soifer, B. T, et al. 2006, ApJ, 638, LI 
Vanden Berk, D. E., Richards, G. T., Bauer, A., et al. 2001, AJ, 122, 549 
Vemet, J., Fosbury, R. A. E., Villar-Marti'n, M., et al. 2001, A&A, 366, 7 
Wang, R., Carilh, C. L., Beelen, A., et al. 2007, AJ, 134, 617 
Weedman, D., Polletta, M., Lonsdale, C. J., et al. 2006, ApJ, 653, 101 
Yan, L., Sajina, A.. Fadda, D., et al. 2007, ApJ, 658, 778 
Yun, M. S., Reddy, N. A.. & Condon, J. J. 2001, ApJ, 554, 803 
Zakamska, N. L., Strauss, M. A., Krolik, J. H., et al. 2006, AJ, 132, 1496 



